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Properties of systems driven by white non-Gaussian noises can be very different from these of 
systems driven by the white Gaussian noise. We investigate stationary probability densities for 
systems driven by a-stable Levy type noises, which provide natural extension to the Gaussian noise 
having however a new property mainly a possibility of being asymmetric. Stationary probability 
densities are examined for a particle moving in parabolic, quartic and in generic double well potential 
models subjected to the action of a-stable noises. Relevant solutions are constructed by methods 
of stochastic dynamics. In situations where analytical results are known they are compared with 
numerical results. Furthermore, the problem of estimation of the parameters of stationary densities 
is investigated. 
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I. INTRODUCTION 

Behavior of many natural systems in contact with their 
surroundings can be described within a stochastic picture 
based on Langevin equations. The basic equation of this 
type reads 



x(t) = f(x)+C(t), 



(1) 



where f(x) is the deterministic "force" representing the 
internal dynamics of the system and ((t) is the "noise" 
describing its interaction with its complex surrounding 
(heat bath). In many cases this noise can be consid- 
ered as white and Gaussian, giving rise to the classical 
Langevin approach used in the analysis of Brownian mo- 
tion. The whiteness of the noise (lack of temporal corre- 
lations) corresponds to the existence of time-scale sepa- 
ration between the dynamics of a relevant variable of in- 
terest x(t) and the typical timescale of the noise. White 
noise can be thus considered as a standard stochastic pro- 
cess that describes in the simplest fashion the effects of 
"fast" surrounding. On the other hand, the Gaussian na- 
ture of the noise is usually guaranteed by assuming the 
surrounding bath being composed of many practically in- 
dependent subsystems and by the fact that the interac- 
tion of x with each of these subsystems is bounded. The 
first assumption allows for considering the noise as being 
a sum of many independent random contributions (in 
thcrmodynamical limit - infinitely many), which math- 
ematically corresponds to the statement that its proba- 
bility distribution is infinitely divisible and stable. The 
second assumption chooses the Gaussian distribution as 
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the only one possessing finite dispersion. However, the 
assumption that the perturbations in the system's dy- 
namics due to interactions with bath are described by 
white Gaussian noise is not always appropriate when de- 
scribing real processes where each of the assumptions 
concerning the noise (e.g. its whiteness or its Gaussian 
distribution) can be violated. In various phenomena in 
physics, chemistry or biology [H, the noise can be still 
interpreted as white (i.e., with stationary, independent 
increments) and distributed according to a stable and 
infinitely divisible law, however, the distribution of the 
noise variable £ is registered as following not a Gaus- 
sian, but rather a more general, Levy probability distri- 
bution. Such situations were addressed for example in 

Refs. SilSiSS&tiltiiEEQEtil- The 

present work discusses some further properties of Levy 
flights in external potentials with a focus on astonish- 
ing aspects of noise-induced bifurcations and explores in 
more detail features of stationary states in Langevin sys- 
tems under the influence of asymmetric Levy noises. 

Levy distributions L aj /3(C; er, /i) correspond to a 4- 
parametrical family of the probability density func- 
tions characterized by their Fourier-transforms (char- 
acteristic functions of the distributions) <j){k) = 
JZo e tkC L a AC, <J, MR being 0, M d 



(j)(k) = exp ikfi — cr Q |A:| Q ^1 — i/3sgn(fc) tan 
for a G (0, 1) U (1,2] and 



7ra\ 
~2~J 



4>{k) = exp 



ik/i — a\k\ ( 1 + i(3— sgn(fc) In \k\ 



(2) 



(3) 
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for a = 1 [l?], EH, 13 ■ Here the parameter a (where 
a G (0,2]) is the stability index of the distribution de- 
scribing (for a < 2) its asymptotic "fat" tail character- 
istics yielding L a ^(<^;a,fi) ~ |^| _ ( 1 + Q ) for large £. The 



2 



parameter a characterizes a scale and j3 G [—1,1] de- 
fines a skewness (asymmetry) of the distribution, whereas 
/i denotes the location parameter. As it is clear from 
Eq. ([2]), Gaussian distribution corresponds to a special 
case of a Levy law for a = 2, with fx interpreted now 
as a mean and a as the dispersion of the distribution. 
However, this special case is somewhat degenerated: the 
dependence of the distribution on j3 disappears due to 
the fact that tan7r = 0, so that all Gaussian distribu- 
tions are symmetric! In general cases strongly asymmet- 
ric distributions (up to extreme, one-sided ones) may ap- 
pear, see e.g. 

El, SI Hi. Such realms are discussed 
much less extensively than the cases of symmetric noises 

BEE El El. 

In many situations, one is interested not in the in- 
dividual properties of the trajectories x(t) but, instead, 
in the one-point probability distributions defined on en- 
semble of trajectories: P(x,t) = (S(x — x(t)))Q at some 
given time t. For stochastic system described by the 
Langcvin equation with an additive white Levy noise 
forcing, the distribution function P(x, t) of the variable 
x(t) fulfills the associated fractional differential Fokker- 
Planck equation. Stationary states (if they exist) can be 
then read from the asymptotic time independent solu- 
tions P(x) = limj-™Ki, t). Such solutions were dis- 
cussed e.g. in [3, yjj [Oa, LLJ] where analysis of Levy 
flights in harmonic/super harmonic potentials has been 
presented. Nevertheless, the discussion presented there 
is far from being complete, since only symmetric Levy 
distributions have been considered. In the present work 
we pay special attention to asymmetric ones. 

In this article we focus on stationary states for a par- 
ticle moving in quadratic, quartic and double well po- 
tentials subjected to a-stable white noises. Theoretical 
descriptions of such systems is based on the Langevin 
equation and/or Fokkcr-Planck equation, which in gen- 
eral is of the fractional order [22j . The model under dis- 
cussion is presented in Section HU Section IIIII discusses 
obtained results, which are divided into three subsections 
regarding results for parabolic and quartic potential (Sec- 
tions IIII Al and IIII B|) and double well potential model 
( Section IIII C[) . The paper is closed with concluding re- 
marks ( Section HV[) . Additional information, regarding 
problem of the dimensionality of the Langevin equation 
is included in Appendix [A] 



II. MODEL 

Let us consider a motion of an overdamped particle in 
a field of a potential force, so that Eq. |T]) takes the form 
of 



Eq. ([4]) can be calculated as [H, 



x(t) = -V'(x) + 0), 



(4) 



a nd C(t) denotes a Levy stable white noise process [1, 0, 
& HH Hj] ■ The value of the stochastic process defined by 



t t 
x(t) = x(0)- I V{x{s))ds+ I C(s)ds 
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x(0) - J V(x{s))ds + J dL a3 . (5) 



Here, the integral f Q £(sWs = /»* dL n jj defines a gener- 
alized Wiener process 0, 0, [23[ that is driven by a 
Levy stable noise, whose increments are distributed ac- 
cording to a stable density with the index a. The Levy 
noise is a formal time derivative of the generalized Wiener 
process. For the time step of integration At, the incre- 
ments of the generalized Wiener process arc distributed 
according to the distribution L a ^(Ax; a(At) 1 ^ a , fi = 0) 
EM [ijj [H, H3] ■ We discuss the overall range of param- 
eters a G (0, 2]; (3 G [—1,1] excluding the case of a = 1 
with (3 7^ 0, for which the numerical results are unreliable 
due to well known numerical instabilities @, [ljj H|| . 
Putting the location parameter of the distribution to zero 
does not influence the generality of our results: Taking lo- 
cation parameter fi to be nonzero is equivalent to adding 
a linear term to the potential (constant drift). Sample 
a-stable probability densities are presented in Fig. [TJ 

The Langevin equation describes evolution of a sin- 
gle realization of the stochastic process {x{t)}. Random 
numbers distributed according to a canonical form of 
characteristics functions given by Eqs. ©-El can be gen- 
erated using the Janicki-Weron algorithm |28l . [29| . More 
details on numerical scheme of integration of stochastic 
differential equations with respect to a-stable noises can 
be found elsewhere @, ©, M, M, HI . 

For a / 1, equation ^ is associated with the follo wing 
fractional Fokker-Planck equation (FFPE) [ISIMIIIIII 
HI 



dP{x,t) 
di 



d_ 

dx 



[fx-V'(x,t)]P(x,t) 



(G) 



d° 



d\x\ a 



P(x,t) 



(i tan ■ 



d d° 



2 dxd\x[ 



-P(x,t) 



where the fractional (Riesz-Weyl) derivative can be un- 
derstood in the sense of the Fourier transform 0, [l2|, 
El M J=rsf(x) = -J-oo ^e~ lkx \k\ a f{k). The frac- 

originate from the form of the 



d\x\° _ 

tional derivatives in 
characteristic function, see Eqs. © and G]), of Levy sta- 
ble variables [l|, [3(J HI, HH. The nonzero asymmetry 
leads to an additional, asymmetric diffusion term includ- 
ing an even, reflection-invariant Riesz-Weyl operator and 
an odd first derivative which changes its sign under the 
x — > — x transformation. The overall order of derivatives 
in the diffusion terms is the same, namely a. 

In the following, the value of the location parameter, fi, 
is set to zero, which guarantees that for a ^ 1 Levy noise 
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FIG. 1: Sample a-stable probability density functions (PDF) 
with a = 1.5 (left panel) and a = 0.9 (right panel). For /3 = 
distributions are symmetric, while for (3 = ±1 they are asym- 
metric functions. The support of PDFs for the fully asym- 
metric cases with fi = ±1 and a < 1 (right panel) assumes 
only negative values for j3 = — 1 and only positive values for 
/0 = 1. Note the differences in the positions of the maxima 
for a < 1 and ot > 1. 



V(x) — > +oo for |x| — > +oo, even the piecewise linear 
one, is sufficient to produce bounded states, i.e., the ones 
with a finite dispersion of the particle's position. On the 
contrary, for the Levy noises with a < 2, the potential 
which grows faster than quadratically in x is required to 
produce bounded states. Furthermore, in the Gaussian 
case stationary probability distributions for a single-well 
potential are unimodal, which is no always true for a Levy 
stable noise with the stability index a < 2 [l2|, Q3 ■ Qual- 
itative and quantitative differences are caused by the fact 
that stable distributions are heavy-tailed and allow larger 
noise pulses with a higher probability than the Gaussian 
distribution [27]|. Moreover, stationary probability distri- 
butions for the additive Levy noises (a < 2), if they exist, 
are not of the Boltzmann-Gibbs type [13, fl3l. [2p|. l37l. |38| . 

In the following sections properties of stationary proba- 
bility distributions for systems perturbed by the general 
Levy noises are discussed. The performed simulations 
corroborate earlier theoretical findings @, [H, [HI, EH- 
Furthermore, the influence of the nonzero asymmetry pa- 
rameter (3 on the shape of stationary distributions is dis- 
cussed. 
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FIG. 2: Exemplary shapes of potential wells used in the study. 
The examination of stationary states of the system perturbed 
by a-stable Levy type noises has been performed for a generic 
double well potential model (V(x) = a; 4 /4 — x 2 /2), as well 
as for parabolic (V(x) = x 2 /2) and quartic (V(x) = a; 4 /4) 
potentials. As discussed in the text, the confinement of tra- 
jectories (observation of bounded states) for a < 2 is possible 
only if the potential slopes are steeper than for a harmonic 
case. 



present in Eq. ^ is strictly stable and standard numeri- 
cal methods of integration of stochastic differential equa- 
tions with respect to a-stable noises (namely the gener- 
alization of the Euler scheme) can be used [2, [l9l |26| . 

The behavior of a system where a particle is subject to 
the additive, white, strongly non-Gaussian noise could be 
very different from the behavior in the Gaussian regime 
[l6| . In the Gaussian case, any potential well such that 



III. STATIONARY STATES FOR A 
"LEVY-BROWNIAN" PARTICLE 

The stationary probability densities P(x) can be ob- 
tained either by analytically solving Eq. ([5]) (which is 
unfortunately possible only for a quite restricted set of 
special cases) or, otherwise, numerically. In such a case, 
there are two approaches possible: either using the dis- 
cretization of Eq. © 0, HH, E3] or employing a 
Monte-Carlo method based on simulation of Eq. ([?]), 
0, 0, US HI]- For the Gaussian noise the solutions 
of the Fokker-Planck equation can be readily obtained 
by using shooting methods and discretization techniques 
[41j . For the general Levy case such solutions can be 
constructed by discretization of Eq. ©, which converts 
partial differential equation to a discrete Markov chain 
|40l. |42|. However, this approach has a drawback of slow 
convergence and possible instability for a < 2 (43[, and 
consequently was not used here. Thus, our data are based 
on Monte-Carlo simulations; our method of solution of 
the Langevin equation is based on the slightly modified 
standard integration scheme for stochastic equations of 
type dH) driven by a-stable Levy type noises [H, [Til [2(| , 
see below. Stationary PDFs were extracted from ensem- 
bles of, typically, N = 10 6 trajectories of a given length 
T max = 10. The value of T max was chosen by trial and 
comparison of numerical estimates of P(x, T max ) for var- 
ious T max (sufficiently long times T max are requited to let 
P(x,T max ) reach stationarity) . A problem related to the 
choice of the simulation time T max is the choice of the 
time step of integration At. The simulations have been 
performed with the time step of integration At = 10~ 3 . 
Such a choice of At guarantees a compromise between 
accuracy and the computational cost of simulations. It 
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is also suggested by earlier studies @, 0, H, [HI, [3 ■ Fur- 
thermore, At = 1CP 3 makes the x-domain in which the 
generalized Euler scheme can be used sufficiently large, 
see Section UlIBl 

For a = 2 (and an arbitrary skewness parameter /3), 
the random force term in the Langevin equation ([4]) rep- 
resents a Gaussian white noise (CMC( s ))a=2 = 28(t — s) 
and the associated Smoluchowski-Fokker-Planck equa- 
tion governing evolution of the probability density P{x 1 i) 
reads 



dP(x,t) _ d 



dt 



Ox 



V'(x)P(x,t) 



dx 2 



P(x,t), (7) 



with a stationary solution assuming the standard Boltz- 
mann-Gibbs form 



P(x) = Afexp 



V{x) 



(8) 



of finite mean and variance. In contrast, Levy flights in 
external potentials exhibit unexpected properties [H, [l2l . 
0, and their stationary PDFs can be shown to pos- 
sess finite first and second moments only, if the imposed 
deterministic forces are derived from steeper than the 
parabolic potentials, see Fig. [I] As an example, in Fig. [3] 
stationary PDFs for a particle moving in the quartic po- 
tential subject to the white Gaussian noise (left panel) 
and white Cauchy noise (right panel) arc compared. Nu- 
merical results were constructed from the ensemble of 
final positions reached after the long time T max obtained 
from the simulation of the Langevin equation (f5]) with 
a = 2, a = 1 and a 2 — 1. cf. Section IlIIBI In forthcom- 
ing Sections we are addressing this point by investigat- 
ing Langevin dynamics driven by asymmetric Levy white 
noises. 



A. Parabolic potential and algorithm testing 

For cr=l, the fractional Fokker-Planck equation ([7]) 
can be rewritten in the Fourier space in the form of 



dP(k,t) 



VP(k,t) - \k\ a P(k,t), 



(9) 



where V is the operator giving the Fourier representation 
of the potential, which can be found in a closed form only 
in the simplest cases, e.g. for polynomial potentials. In 
the case of asymmetric distribution the analogous equa- 
tion reads 



dP(k,t) 



dt 



VP(k,t)-\k\ a 1 - i/?sign(fc) tan ■ 



P(k,t). 
(10) 

The choice of the parabolic potential V(x) = x 2 /2, see 
Fig. [TJ results in its Fourier transform V = —k-g£. For 
symmetric a-stable noises, the corresponding equation 
for the stationary PDFs is QJ, d 



clP(fc) 
dk 



-sign/c|fc| Q - 1 P(fc), 



(11) 





FIG. 3: Stationary solutions to the Smoluchowski-Fokker- 
Planck and to the fractional Fokker-Planck equations associ- 
ated with the diffusion (a = 2, left panel) and superdiffusion 
(a — 1, right panel) in a quartic potential V(x) = x 4 /4. For 
a = 2 numerical results were constructed by discretization 
techniques and shooting methods [4l]]. Simulation parame- 
ters: N = 10 6 , At = 10" 3 , iVbms = 100 and T max = 10. 



and its solution in the Fourier space reads 
P(k) — cxp 



(12) 



i.e., the stationary solution is a symmetric Levy distri- 
bution, see Eqs. © and ([3]). Consequently, the vari- 
ance of the stationary solution is infinite. Therefore, the 
parabolic potential is not sufficient to produce bounded 
states for a particle subject to the action of a Levy 
noise [1 [12, El, 0, ISl. For potentials steeper than the 
parabolic well, the confinement of superdiffusive trajec- 
tories becomes possible but, additionally, the additive 
Levy white noise could induce bimodality in the station- 
ary PDF. 

For general a-stable driving (a ^ 1) the stationary 
solutions obey the equation 

= -signfc|fc| Q - 1 P(fc) +z/3tan — Ifcl^P^), 
ok 2 

(13) 
(14) 



and its solutions is 
P(k) = cxp 



k\ a / . . 7ra\ 
\ \ — zpsign(fc) tan — 



2 / 



For the nonzero asymmetry parameter, /3, like for the 
symmetric noise, the stationary probability density func- 
tion is a stable law with the same stability index a and 
the asymmetry parameter j3 and a different scale param- 
eter a' = aoT 1 ! - (here a' = q _1//q ). The location pa- 
rameter fx of the resulting distribution is zero. The ex- 
istence of these analytical results allows to use the case 
of parabolic potentials as a test bench for our simulation 
algorithms. 
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Thus, for the testing' purposes, large samples of long 
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are observed for the location parameter \x. In some situa- 
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were compared. Estimated values of parameters for both 
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values of T max are consistent. Therefore, only results for 
Tmax = 15 arc presented, see Tabs. HI HTT1 

TABLE II: Continuation of Tabled (a = 1-1)- 
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TABLE I: Theoretical and estimated values of stationary PDF 
parameters. In every cell numbers in bold indicate theoretical 
values of parameters, followings rows estimated parameters 
using quantile evaluation (2 nd row), characteristic function 
evaluation (3 ld row) and maximum likelihood method (4 t h 
row). Distributions' parameters were estimated by use of [441 ] . 
Simulations' parameters: T max = 15, At = 10 -3 . Samples 
contain not less than 10 6 elements. * indicates cases when the 
software applied warned about some problems in estimation 
of sample parameters. 



TABLE III: Continuation of Table|T](a = 1.8). 



B. Quartic potential 



For the quartic potential V(x) = x 4 /4 (V = kJ^j) 
and symmetric a-stable noises, fractional Fokker-Planck 
equation in the Fourier space has the form 



d 3 P(k) 
dk 3 



sign/fc|fc| Q ~ 1 P(fc). 



(15) 
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The solution of Eq. (|T5|) is known for a = 1 [12], [L 



v/3 



exp 



/ V3|fc| 



, and the cor- 



responding stationary solution in the real space reads 



Pa=i(a:) = 



1 



7r(l - a; 2 + a; 4 ) ' 



(16) 



The stationary solutions (fT6|) of Eq. (p~5|) have two main 
properties. First of all the stationary PDFs are not of 
the Boltzmann type, as typical for the stationary states 
for system driven by Levy white stable noises with the 
stability index a < 2 [2(j. Additionally, the stationary 
probability density function for a quartic Cauchy oscil- 
lator is bimodal, with extremes located at x = dsz 1 /^ 
ElEE El- Fig.[3]prescnts stationary PDFs obtained by 
the simulation of the Langcvin equation along with 
theoretical lines for the Gaussian (left panel) and Cauchy 
(right panel) quartic oscillators. Moreover, the parabolic 
addition to the quartic potential V(x) = ax 2 /2 + x 4 /4 
(a > 0), can diminish or even destroy the bimodality of 
the stationary PDF [3 El 03 ■ Finally, for the general 
a-stable driving and the quartic potential the stationary 
density fulfills 



d 3 P(k) 
dk 3 



signfcl/cl"- 1 P(k) 



i/3tan^|fc| Q - 1 P(fc). 



(17) 

Quartic potentials pose some additional difficulties to 
Monte-Carlo simulations making it necessary to slightly 
modify the standard techniques of integration of stochas- 
tic differential equations driven by a-stable Levy type 
noise • Due to heavy tails of stable distribution large 
random pulses are much more likely to occur than in the 
Gaussian distribution leading to very long jumps from 
time to time putting a particle to the position where the 
force acting on it is very large. In this case approximat- 
ing the deterministic drift by — AtV'(x) is too inaccu- 
rate whatever small time step is chosen. It can result in 
switch of the particle to the other side of the origin in 
such a way that new particle's position is more distant 
from the origin than the initial one leading to numerical 
instability and escape of the particle to infinity Of 
course, this effect is weaker when smaller integration step 
is used. However, taking smaller steps was proven not to 
give an effective solution to the problem. Our approach 
to it is based on separating noise and deterministic drift 
and integrating the last one analytically, by solving the 
differential equation x(t) = — V'(x) to obtain x(t + At) 
for a given initial condition x(t). Such a step (involv- 
ing the solution of an algebraic equation) is more time- 
consuming than the Euler integration step and is abso- 
lutely superfluous for small and moderate x. Therefore, 
exact integration of the deterministic part is performed 
only for large x, \x\ > x tr , while the noisy part is always 
integrated in the standard Euler way. In the simulations 
we took xtr = 15 as motivated by analytical estimates 
and by numerical tests. For testing purposes constructed 
numerical results for the Cauchy noise (a — 1) were com- 
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FIG. 4: Influence of the time step of integration and duration 
of simulation on portion of non escaped trajectories at time 
Tmax for V(x) = a- 4 /4 — x 2 /2. x tr represents threshold value 
of x such that for |x| > str deterministic part of Eq. Q 
is integrated analytically. The case x tT = 10 12 corresponds 
simply to the standard Euler scheme. For more details see 
the text. 




FIG. 5: Stationary states for the generic double well potential 
model subjected to the a-stable driving with a — 0.5 and 
various f3. Simulation parameters: N = 10 , At = 10 -3 , 
M>ins = 100, xtr = 15 and T ma x = 10. Stationary states for 
totally skewed noise, i.e., = ±1 are one-sided. 



pared with the known analytical solution, see right panel 
of Fig. [3l leading to the excellent level of agreement. 



C. Double well potential 

The results for the double-well potential model, see 
Fig. [1] were constructed by the numerical method de- 
scribed in the previous section (Section IIII Bp . Further- 
more, we compared the influence of decreasing time step 
of integration and x tr . The results of comparison of both 
methods of reduction of number of numerical escapes are 
summarized in Fig. 3] where the influence of the time step 



FIG. 6: The same as in Fig. [5] for a = 1.1. 




FIG. 7: The same as in Fig. [5] for a = 1. 



of integration At, duration of the simulation T max on the 
portion of valid (non escaped) trajectories are compared. 

Stationary solutions shown in Figs. [SH7] are obtained 
for the generic double well potential model, i.e., V(x) = 
x 4 /4 — x 2 /2. In the simulation the whole allowed range 
of a and (3 was examined, in figures only a limited choice 
of representative values of noise parameters is presented, 
namely the same as ones used in Tabs. U tillll For a < 1 
with \(3\ = 1 stable distributions are one-sided, there- 
fore, stationary solutions for a = 0.5 with (3 = ±1 are 
different from zero only on the one side of the origin. 
For \(3\ < 1 stable distributions takes all real values as 
manifested by nonzero probability for all x, see Fig. [5] 
Furthermore, the symmetry of noise and the potential 
is reflected in the symmetry of stationary densities, i.e., 
solutions for —j3 can be constructed by the reflection of 
the solutions for /3, see Figs. [SHU For (3 = with any 
a stationary densities are bimodal and symmetric along 
x = 0. For a > 1, the support of stable densities as well 
as the stationary distributions is whole real line. Conse- 
quently, even extreme values of the asymmetry param- 
eter (3 = ±1 are not sufficient to switch the probability 



FIG. 8: Median (quantile go. 5) of the stationary probability 
density P(x) (upper panel) and asymmetry of stationary dis- 
tributions measured as the fraction of the probability mass 
on the left hand side of the origin, i.e., P(left) = P(x)dx 
(lower panel). For symmetric noises resulting stationary den- 
sities are symmetric along x = and consequently P(left) = 
1/2 and qo.5 = 0. Different symbols correspond to different 
values of the stability index a: 'o' a = 0.5, 'x' a = 1.1 and 
'□' a = 1.8. 



mass to the one side of the origin, see Figs. |5] [3 Further- 
more, it is well documented in the lower panel of Fig. [8] 
where P(left) = f®^ P(x)dx is presented. Theoretical 
considerations as well as the probability of being in the 
left state, P(left), indicate that stationary densities can 
be one-sided only for some totally skewed a-stable noises 
with small a, i.e., (3 = ±1 with a <C 1. For a > 1 with 
(3 = ±1 two maxima of stationary PDFs are visible. In 
upper panel of Fig.[8]locations of the median value, which 
may be considered as the next measure of asymmetry of 
stationary probability densities, are presented. 

Very small values of a (a < 0.5) pose special difficulties 
for simulations. Our simulation of Eq. (j4|) starts with the 
initial condition x(0) = 0. The initial transient peak of 
the probability density at this value is rather persistent 
for small a, so that the simulation time has to be long. 
On the other hand in this case simulations are prone 
to escape of trajectories to "infinity" due to too strong 
noise pulses and require very small At, so that the overall 
quality of such results is not very good. Therefore, the 
results for a < 0.5 are not presented here. 

Another method to present the results for stationary 
distributions is the use of the effective potential. In gen- 
eral, the same stationary probability densities that are 
recorded for the double well system driven by a-stable 
noise, see Eq. can be observed in the Gaussian regime 
in the effective potential V c s(x) = — lnP(x). Sample ef- 
fective potentials corresponding to stationary states for 
a = 1.1 from Fig. [6] are depicted in Fig. [9l The same 
stationary solutions can be observed for motion in a sim- 
ple potential, like double well potentials, and a-stable 
stochastic driving or for potentials of the complicated 
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FIG. 9: Effective potentials V c g(x) = — \nP(x) correspond- 
ing to stationary densities from Fig. [6] The effective poten- 
tial V a s{x) together with the white Gaussian noise results in 
stationary densities of the model Q which are presented in 
Fig.[6] 



form, see Fig. [9l and standard white Gaussian driving. 
Despite the fact that stationary states are the same, other 
characteristics of these two processes arc different. 

IV. SUMMARY AND CONCLUSIONS 



potential: the asymmetric densities correspond to asym- 
metric potentials, i.e., to deterministic forces which break 
the symmetry of the initial problem. The Gaussian noise 
itself is always symmetric. For systems driven by Levy 
noises (a ^ 2), an asymmetric stable noise together with 
symmetric static potential is sufficient to produce asym- 
metric stationary densities. In this situation the asym- 
metry of stationary states originates from the asymmetry 
of the stochastic driving and can be controlled by chang- 
ing the parameters of the noise. 

Our main studies have been performed for the generic 
symmetric double well potential model. Here the asym- 
metry of the stationary distribution (as measured by the 
probability of being in the left/right state or location of 
the median) was investigated as a function of the param- 
eters of the noise. The asymmetry of stationary state 
decreases with increasing a, see Fig. [8] Finally for a = 2 
with any value of the asymmetry parameter (3 Gaussian 
scenario is recovered and stationary density is fully sym- 
metric. We also checked whether stable asymmetric noise 
can produce unimodal stationary PDFs in the double well 
potential. Such situation can indeed be observed for fully 
asymmetric noises = 1) with a <C 1, e.g. for the 
Levy-Smirnoff noise (a = 0.5, (3 = 1). 

APPENDIX A: DIMENSIONALITY OF THE 
LANGEVIN EQUATION 



In the present work we investigated the form of sta- 
tionary probability densities of a position of a particle 
subject to a deterministic potential force and to a Levy 
noise, paying special attention to the case of asymmet- 
ric stable noises. Stationary density functions for system 
driven by a-stable Levy noises (a ^ 2), if they exist, are 
not of the Boltzmann-Gibbs form. 

For parabolic potential the stationary density func- 
tions are a-stable laws with the same stability index a 
and asymmetry parameter (3 as ones of the noise. The 
only difference is the scale parameter of the resulting dis- 
tribution. Therefore, this case can serve as a benchmark 
for our simulation algorithms. By use of the special soft- 
ware [H| for estimation of stable law parameters, the 
parameters of stationary densities have been evaluated 
leading to a very good level of agreement between theo- 
retical and estimated values of parameters. Natural con- 
sequence of the Levy type of the stationary densities for 
parabolic potential is divergence of the variance of the 
particle position. To produce bounded states, i.e., states 
with finite variance of a position, potentials steeper than 
quadratic are necessary. Thus, for the quartic potential 
variance of stationary densities is finite. Furthermore, for 
white Cauchy noise analytical solutions to the stationary 
problem is known. Here again, numerical results fully 
agree with earlier [l2|, (HI theoretical findings. 

In the case of Gaussian noise (a = 2) the symmetry of 
stationary density (being the Boltzmann-Gibbs equilib- 
rium distribution) reflects symmetries of the underlying 



We have investigated the dynamic stochastic process 
modeled by the (overdamped) Langevin equation of the 
form 

x(t) = — +D 1/a ((t) = ^+a((t) 7 (Al) 

where: x - is a position of the particle, 7 - stands for a 
friction coefficient, m - is particle's mass, D (er) - repre- 
sents strength of the noise and £(i) - is Levy stable white 
noise characterized by the stability index a (a & (0,2]) 
and the asymmetry parameter (3 ((3 6 [—1, 1]). The force 
acting on a particle is determined by the external poten- 
tial, f(x) = —dV(x)/dx. 

Corresponding units in Eq. (|A1[) are: [x] = [length], 
[7] = l/[t], lf(x)} = [V'(x)] = [m] x [length]/^] 2 = 
[force], [V(x)] = [force] x [length] = [energy], [D] = 
[length] a /[t] ([a] = [length]/^] 1 /") and [£(*)] = 
~ 1//q . Stability index a and asymmetry parame- 
ter (3 are dimensionless. In the asymptotic limit of a = 2 
the Levy white noise is equivalent to the Gaussian white 
noise and it has standards units, i.e. [Ca=2(i)] = [£(£)] = 

VVWl 

By the set of transformation: t — > t/to and x — > x/xq. 
Eq. (|A1[) can be transformed to the dimensionless form 

x{t) = f(x) + D 1 ^^) = /(*) + <r<V), (A2) 

which is of the same type like Eq. ^ because the 
(rescalcd) noise intensity can be incorporated to the dis- 
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tribution of the particle's position increments. An alter- 
native way of getting dimensionless form of Eq. (|A1|) can 
be found in a recent work [47j |. 

For the single- well potential V(x) = ax n /n 



x 



(A3) 



Thus V(x) — > x n /n and a — ► 1. Therefore, for the sin- 
gle minima potential, the only relevant parameters are a 
(stability index) and j3 (asymmetry parameter). 

For the generic double-well potential V(x) = —ax 2 /2+ 
bx 4 /4: 



a 



7m 



(A4) 



Thus V(x) 



-x 2 /2 + x 4 /4 and a — » at^ a /x . In con 



sequence, for the double-well potential, the only relevant 



parameters in a dimensionless form of Eq. (|A1[) are the 
Levy noise parameters a and j3, and a (rescaled) noise 
strength. From the above analysis, it is clear that the 
choice of xq , io introduce scales to the system which are 
directly related to its dynamical parameters. 
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